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ABSTRACT: The Feynman- Vernon influence functional formalism is used 
to derive the effect of a black body, treated as an environment, on a massless 
scalar field in 1+1 dimensions. The black body is modeled as a finite region of 
space in which an independent ohmic heat bath is weakly coupled to the field 
at every point. A weak coupling approximation is developed which implements 
the concept of an ideal black body in the context of quantum field theory. The 
calculation takes advantage of the suitability to harmonic oscillators and free fields 
of Bargmann-Fock coherent state variables, whose convenience is illustrated by a 
preliminary derivation of the thermalization of a single harmonic oscillator by a 
heat bath with slowly varying spectral density. The black body model exhibits 
absorption, thermal equilibrium, and emission consistent with classical results for 
black bodies. It is argued that this model is in fact a realistic description of a 
very fine, granular medium, such as lampblack. 
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phenomenon will be shown in this paper to explain the fact that the thermal 
radiation of a black body, although it is generated via the coupling of photons 
to charged matter, is not proportional in intensity to the weak electromagnetic 
coupling^ . 

Second order environmental effects of these kinds are not described by clas- 
sical external sources, but may be accurately expressed in the influence functional 
formalism of Feynman and Vernon[9]. This formalism has already bccu used to 
examine quantum measurement [6]. In the present paper, it is applied to radiation 
and absorption by an ideal black body of finite extent. Inasmuch as we are able 
to describe radiation propagating away from a heat source, our results go beyond 
the standard equilibrium treatments. 

This paper is organized as follows. The problem of a simple harmonic os- 
cillator driven to thermal equilibrium by contact (velocity coupling) with a heat 
bath is presented first, in the comparatively little-used formalism of the path inte- 
gral in Bargmann-Fock coherent state variables. Although these variables arc less 
familiar than positions and momenta, they are actually particularly well-suited 
to problems involving harmonic oscillators, including any calculations involving 
free quantum fields. They are introduced here in a demonstration that standard 
thermal equilibrium is established, after a sufficiently long time, whenever the 
coupling to the bath is weak enough and its spectral density varies slowly enough; 
a spectral density obeying a particular "ohmic" power law is not required. For 
simplicity, we then assume an ohmic environment with constant spectral density, 
and consider a massless quantum field coupled to a heat bath restricted to a finite 
region of space. Here the coherent state variables are particularly advantageous, 
allowing us to clearly see that this model exhibits uniform absorption and black 
body radiation. We then conclude with a discussion of the realism of the model 
chosen, the issues raised by our analysis, and some related problems for future 
study. An appendix reviews the important justification, originally given by Feyn- 
man and Vernon [9], for treating a generic weakly-coupled environment as a bath 
of independent harmonic oscillators. 



It is a fact that should be more surprising than it is, that the Stefan-Boltzmann 
constant a = j^g^aj is independent of the elementary charge e. 



II. OSCILLATOR WITH 
IN BARGMANN-FOC 



A. The Bargmann-Foc 

When using path integrals to obtain 
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Consider the following coherent states as basis elements for bras and kets: 



where the \n) are harmonic oscillator energy eigenstates, and a and a are complex 
numbers (not necessarily complex conjugates of each other). Since a whole 2- 
planc of parameter space is larger than needed to span the states of an oscillator, 
we constrain a to lie on a line C through the origin, and a to he on another 
line C through the origin, such that C is perpendicular to C. We then have an 
expression for the identity operator, 

^ [ da [ dae-""|a>(a| =1. (2) 
2^' Jc Jc 

(This relation may be verified by inserting the definitions (1) and using differen- 
tiation under the integral sign, treating the expression as a distribution.) 

Since \a) is an eigenstate of the annihilation operator a, any operator ex- 
pressed as a normal-ordered function of creation and annihilation operators has 
a simple matrix form in this representation: 



(a| ■.A{d\a): \a) =A{a,a) . 



(3) 



In particular, for a Hamiltonian H =:H{a\a):, we have the following matrix 
elements of the infinitesimal time evolution operator: 



(a|e-'"-^|a) =e- 



i5t /i(a,a) 



+ 0(5^) 



(4) 



(Throughout this paper wc set Ti = c = 1.) Using (2) and (4), we can derive the 
path integral for a transition amplitude in the usual way, by inserting identity 
operators between an infinite succession of infinitesimal time intervals. We obtain 
the Bargmann-Fock path integral 



{af\e 



i f ds 
Jo 



H(s), 



|Qi) = / VaVa5{a{0) - a,)5{a{t) - a,) 



X exp 



afa{t) + a{0)ai 



+ i 



f 

Jo 



ds 



2i 



H{a,a] 



(5) 

s being used as a time parameter, with a = ^q(s). Note that, according to the 
delta functions in (5), a{s) and a{s) have only one boundary condition each, at 
the initial time s = 0, and the final time s = t, respectively. 



In the case where the Hamiltonian H is 
up to a normalization factor, to the value of 
extremized. As usual, this condition is satisf 
equations. In this case, the Euler-Lagrange e 
that only the boundary terms do not vanish 
therefore simplifies to 

{af\e Jo \a,) = — e^v * 

for some normalization factor Z. The subscri 
endpoint values of the solutions of the Eulei 
in the case of a harmonic oscillator driven bj 

where / and / arc c-number functions reprt 
these solutions may be obtained by direct 

Green's functions. The resulting endpoint ' 
the driving sources, as well as on the initial < 



a{t) = OiC -I- 



d 



q(0) = Q/e + i I a 
'o 

While this paper deals with interactions thi 
external sources of this type, similar calcul 
Euler-Lagrange equations, and boundary v; 
below. 

Because {n\a} = -^=, it is easy to obtai: 
states, without integrating, from (6). One sii 



(af\e 



i f ds H(s) I 
Jo 



and extracts the co-cfEcients of the powers of 
oscillator, this gives a clear interpretation of t 
The terms in a{t) proportional to a*, and i: 
energy-conserving evolution of initial states 
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a{t) describes excitation of the system, by the driving source, while the other 
term in a(0) implies absorption by the source. This interpretation is emphasized 
here in order to provide some intuition for the results of Section III, in which a 
quantum mechanical black body will appear instead of a classical driving force, 
but in which the action will still be quadratic. All of the system's dynamics will 
still be expressed in the boundary values of a and a variables. 

The basis in which a path integral is expressed determines the type of bound- 
ary conditions which the paths must satisfy. For given boundary conditions, 
though, one is free to change variables in the path integral. For example, we are 
free to re-express (5) as 

(a/|e Jo '< 



involving free quantum fields. We assume th 



/ 



VRDQ 5{flQ{0) + iP{0) - V2^iai)S(p,Q{t) - iP{t) - 
at[Q.Q{t) + iP{t)] + a4nQ(0) - iP(0)] 



2naf) 



X exp 



(10) 



+ i 



ds 



PQ 



2V2n 

QP 



-h{P,Q) 



-.{QQ-iP) 



using the transformation 

^ -.{nq + iP) 



(11) 

and absorbing the constant Jacobian in the new measure. While the integral now 

involves familiar-looking P and Q variables, the transformation we have used will 
in general make P and Q complex, and the usual P{t) and Q{0) boundary con- 
ditions do not apply. They are replaced by conditions on the linear combinations 
a and a, reflecting the fact that (10) is still a transition amplitude between co- 
herent states. The P, Q type of variables will be more convenient for solving the 
equations of motion of Section III, while a, a variables will still be better suited 
to the initial and final states. In the path integral formalism, there will be no 
reason for us not to take the best of both worlds. 



B. Oscillator coupled to a heat bath 



Consider now a system consisting of a simple harmonic oscillator of natural 
frequency f2, minimally coupled to a heat bath made up of environmental oscil- 
lators. This is a toy model, but it is very relevant to more significant problems 



H=Hs 



HsHo =i(P'-FO"Q') =0 



Henv=-^ J (LjI{uj)\^{j 



where a = -^^{^Q + ^-P); and g is a couplin; 
n, and will use this fact extensively below, 
density, whose properties will be further (lis 
the operators, (12) represents the choice of 
of standard electromagnetism. This choice 
H positive definite and of reducing the size 
normalization of Q. 

We are interested in the behaviour of t 
we must integrate out of the problem all of th 
Assume that the initial density matrix for th 

p(0) = p(0) ® 

where p{0) acts in the Hilbert space of the 
initial density operator for the environment, 
observed oscillator at time t is then given, 
propagator equation 



p{a;,a'f,t) = {aj\p{t)\a'f) 

daida'id£,d^ 



J {af , a'f , a. 



This equation may be translated into the m 
basis, by defining the propagator Ktimn in t 
that 

J{af,a'f,ai,ai;t) = Kk 

klmn 

which implies that 



.(t) =^p.i(0) 



k,l 
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(17) 



J(q/, a^; a'i, at; t) may be obtained explicitly from a path integral similar to 

(10): 

X 5{a'{t) - a'f)5{a'{0) - e*(^+^+^^ . 
In this expression, A is the action term 

A=^ [ dsPQ-QP- P'Q' + Q'P' - {P^ - P'^) - Q\Q^ - Q") ; (18) 
^ Jo 

B is the boundary term 



B ■: 



oifait) + Q;(0)a. + q-q'(O) + a'{t)a'f 



2i 



(19) 



whore the a variables are given in terms of the Q's and P's by (11); and V — 
V[Q,Q'] is the influence phase. (The quantity e*^ is known as the influence 
functional.) 

V contains all the information about the heat bath that is relevant to the 
observed oscillator. It may be computed in a straightforward manner. (This has 
been done using the Lagrangian form of the environmental path integral: the 
Lagrangian obtained from Henv is 



i I I{iu)duj(^ql+2gQq^-u\l^ 



(20) 

which is that of the velocity coupled model analysed in Reference [11].) In the 
case where the initial state of the bath is thermal at inverse temperature /3, one 
finds that 

V[Q,Q'] = 

2 poG 



]9_ 
2 



■ j dujijjl{uj) J '^^ J ds' {Q — Q')s(^{Q — Q')s' coth.^^ cos tij{s — s') 

- i(Q + Q')s' sinw(s - s')^ 

2 /"oo rt 

j dwl{^) J ds{Q^-Q'^)s 
I ds I ds'{Q-Q'),({Q-Q')s'Hs-s')-iiQ + Q'),'vis-s')^ 



-1-2 

''2 



Jo 

2 ft 



dsiQ^-Q'^) 



We use the notation (Q - Q'), = (Q(s) - Q'{ 
define the quantity and the functions r? ai 

The exponent i{A + B + V) in (17) is qu 
constant factor to the extremized value of th( 
subject to the constraints imposed by the bo 
order Euler-Lagrange equations for the fou: 
Q±Q', P±=P±P': 

Q.{s) = P.{s) 
Q+{s) = P+is) 

P.{s) = -{Q' + f/)Q.{s) - 



-{Q' + fi^)Q+{s) + g' 



+ ig^ / ds'Q-{ 



P+{s) = 



These equations together imply that A + V 
duces to one of endpoint values at the extrer 
1 



Jiaf,oi'f,a'i,ai;t) = — 



Z 



^i{A+B+V) 



a/a(t) + a(C 



(21) 



= -exp. 

where Z is a normalization constant. Thus 
(25), subject to the boundary conditions 

i(n(Q++Q-) + i(P+ + J 
^{Q(Q^-Q_)-i{P+-F 
i(n(Q+-Q_)+i(P+-/ 

^{i}{Q+ + Q-)-i{P++f 

imposed by the delta functions in (17), and 
that these solutions determine. 

The easiest way to solve these coupled 
work in Fourier space. The problem then re( 
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of a complex resonant frequency, as a pole in Fourier space. To begin, we dif- 
ferentiate (22) and substitute it into (24), obtaining a decoupled, second order 
integro-differential equation for Q-: 



Q. + (fi' + M^)Q_ = / ds' Q-{s')ri{s - s') 

J s 



(28) 



We then define the Fourier transform Q-{s) = ^ J^^dujQ^e^'^^* ^\ and use 
the definition ri{s — s') = J^dut ujI{u!) sin uj{s — s') (implicit in Eq. (21)) to obtain 



doje 



iu}(t—s) 



-0.^(1 +fl^A(a.O) ) -|-a;/(|u;|)P 



djjj 



LO' — OJ 



(29) 

where we have combined the jj? = duj I{oj) term with the contribution from 
the lower limit of the s' integral to form the frequency re-normalization term 



Jo 



du! 



(30) 



The P in front of the integral sign denotes the Cauchy principal value. We assume 

that \{u)'^) « jj. 

The crux of the problem now appears to be the Hilbcrt transform in the 
rightmost term of (29). The Hilbert transform has the property that 

1 



/ LJ' — Ul UJ' — (x 



X + iy) sgn(j/) ^ _ _|_ 



With this in mind, we choose the ansatz 



QZ = 



(31) 



(32) 



w — (n + 17) 

for some Q, and 7 > to be determined, and A a constant to be fitted to boundary 
conditions. Equation (29) then becomes 

°'"*n~n ('"'^^ ^ ^'^^ " " = ■ (33) 



diu — 

-00 ^^-(n + i7) 



We now invoke our assumption that -f^ is a small quantity, and require 
further that the spectral density /(w) is always much smaller than and is 
slowly varying near u! = Q.. With these conditions, we can choose 

„2 



7=^/(0) 



n = 



(34) 



and take the chosen in (32) as the leadin 
2 . 

The first order correction to Q-{s) will 



f 

J — c 



iw(t-s) 



-(/(0)-/( 

1 ^ 

51 — 17 + — 47/ 



■0-17 V 2 
1 



The poles at ±Q + i-y contribute "correctio: 
zeroth order solution. For a wide range ol 
no sharp peaks near ui = Q, the remaining 
particular, we can have a gradual cut-off at 
realistic spectral densities. 

We can substitute Cl —> —Cl and the p 
With our assumptions about and /(w), the 
to (28) 

Q_(s) ~ Af/'^'-^'e^^""*' + 

valid up to corrections of order ^. 

This solution is anti-damped, growing 
haviour is usually considered unphysical, bi 
such solutions do not apply in this case. De 
similarly to a position variable, it is really a 1 
onality in the position basis of the oscillator 
represents the oscillator's mean displacemen 
classical equation of motion, but rather a to 
quantum mechanical system. In classical me 
conditions must apply at the initial time, tl 
results from a "runaway" solution; but iu qi 
final time, with its own boundary conditior 
be made to conform. And in this case, the a 
essential feature responsible for establishing c 
has long been known to occur in the limit wl 

This series may bo constructed by replacing 
delta functions. (A continuous spectral densi 
tion for such a discrete spectrum, anyway.) 
(ai-aj^)^+e^ then allows us to use (31) repeat 
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We now turn to Eq. (23), differentiating it and applying (25). We obtain 
the damped and driven oscillator equation 

(5+ + (n^ + M^)Q+-/ / ds'Q+{s')r){s-s')=ig^ [ ds' Q-(s>(s - s') ■ (37) 
Jo Jo 
The homogeneous part of (37) is simply the equation formed by substituting 

(t — s) — > s in (28). We therefore have the free (approximate) solutions 

Q° (s) ~ Be-'^'e""' + Be'^'e-"" . (38) 



The source on the RHS of (37) is given by (36) as 



' fds' 
Jo 



f 

J —I 



Q-{s')i/{s — s') 



w coth '^I{\uj\)e'"'' { A- ^ + A 



7)t ^—iujt 

17 



a; — n + 17 



(39) 

For the particular solution, we presume that Q+{s) does behave as a damped 
oscillator, and therefore consider the candidate 



2 r'^ ,.,T(\oj\)coih^e'^' 

Q, — i'y) (w — n — i'y) 



( 



X A- 



oj + f2 + 17 



-iuit (-in — r)t \ 
+ A^ i— ) (40) 



o — c 



di^e'-'^\AQt^ AQt) 



We examine first the term proportional to A in the LHS of (37): 

dwe*"' \Qt{^^ - a;'(l + /A)) - |-w/(|a;|)P / dw' j 

-OO O — OO ^ 



7 *l 

do; e 



Q::^(n'' -w^(1 + /A)) - ^a;/(|a;|)ncoth^7(|n|) 



X P 



duie'' 



io' - Lu (io' + Q + i'y){uj' + n - i7)(a;' - - 17) 
Q^(n^-a;^(l + ff^A)) 



2 /(|a;|)ncoth^(e('^-^'* -e-'"*) 



2 



da;e" 



(w + + 17) (a; + n — 17) {u! — Q — i'y) 
^ g(»n — r)t _ g-iwt 

a; + n + 17 ' 



(41) 



neglecting terms of 0(^5-). In (41) we have U 
sharply peaked around u! = —Cl, and that un 
region. We have also invoked similar appro: 
(36). 

Adding the similar term proportional to . 
we can see that (40) does indeed satisfy (37) 
(40), once again using the fact that 
that 



f 

O — C 



„(iri-7)(t-s) 



■ iTr/n coth ^7(0) { ^ 

" 2 ^ -^V 4z7n 



But the decaying terms in (42) are solutions 
may be absorbed into (38). Including the si 
a change of sign because coth ^ is odd), \ 
solution 



Q+(s) ~ (Be-*^'+Be'^'')e-^'-coth^(A 



Using (36) and (43) in (22) and (23), wi 



iP_(s) ~n(yle*^(*-">e^f"-*' - A 



iP+(s) ~-^coth^(^e'^(*- 



+(Be 



-if2s 



Be*"")e" 



We can now impose the boundary conditions 



2fiai ~ fi((l-coth^) 



2Q,a; ~ n(-(l + coth — 



B^. 



f^ois ~ r2((l + coth 

pKiol, ~ n((coth^ - 1) 
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In the thermal Umit where e ^ and > 0, these imply 



A = 



ficoth^ + l 



B 



Using (46) and (11), we now obtain 

a{i) = e~^^af 
a {t) = e~^^a.f 



a'(0) = Qi 
a(0) = a\ 



Applying these results to (26) , we find the propagator 



(46) 



(47) 



0. Hence the distinction between O 



state, within the limits of our approximation 

We therefore conclude that a harmonic 

pling to an environment with slowly varying t 
to thermal equilibrium on the thermal time i 
of Reference [7] , where the path integral in p 
thcrmalization for a weakly coupled cnvironrr 
our velocity coupled model, this would be eqi 
[7] notes that more general dissipative envir( 
but does not treat them explicitly.) We co 
that, in the case of the single oscillator cou 
scale that makes the thermal effects of the I 
weak coupling and slowly varying spectral i 
environment from communicating to the sys 
temperature, and so guarantee the canonical 



J = exp (a'a; + e ''"a/a'/) 



(48) 



whore Z will turn out to be the partition function of the canonical ensemble. 
Using (14), we find that the final state of the oscillator has the density matrix 



p(a,a) = — exp(e '^^aa) , 



which is the canonical ensemble in coherent state variables, regardless of the initial 
state. We may translate this result into the standard energy representation, using 

(15) to discern 



which implies that 



(t) = (1 - e-'^")5.,e-'"'^"5^„ + OC-) + 0{e-~") , 



p^„(t) ~ (1 - e-''")5^„e-"'''" . (49) 

- 2 

(Although n = 51 + 0(^2-), one might be concerned that for low temperatures, 
where P becomes arbitrarily large, /3(jr^A(n^) would not be small, so that it would 
be wrong to write Q instead of Q, in (49). But because ff^ « Q,, regardless of 
temperature, fig^ can only be non-negligible when is so large that e~^^ ~ 



III. BLACK BODY 1 



A. Setting up th 



We now present our main result. Whei 
to an environment which is confined to a finil 

lengths long compared to the coupling tend 
this region, but shorter wavelength quanta i 
When the thickness of the region is large c< 
sion through it is also negligible, and the q 
as thermal radiation. After the initial time 
ronment are uncorrelated, this radiation prop 
region; but behind a wavcfront region whose 
square of the coupling, it is time-independe 

Cooling of the emitting medium is not consic 
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(i.e., one with slowly varying spectral density /(w)) may be analysed by straight- 
forwardly extending the approach of Section II, for simplicity in this Section we 
consider an ohmic spectral density (/(a;) = 1). The environment wc discuss thus 
represents a black body, absorbing all frequencies (above the reflection regime) 
indiscriminately. 

Also for simplicity, we consider a scalar field in 1 + 1 dimensions; the ex- 
trapolation to a vector field in 3 -I- 1 dimensions is straightforward enough, but 
will only be touched on in this paper. We let our spatial co-ordinate x £ [—A, A], 
and impose periodic boundary conditions on the field 0. (We will eventually let 
A oo, of course, but this infrared regulator will turn out to be convenient.) The 
black body is placed at the origin, filling the region — L < x < L with a uniform 
medium consisting of independent ohmic heat baths at every point, minimally 
coupled to the field*. The Hamiltonian is 

f-A 

dx{ll^(x) + d^<j>\x)) 



«4 



+ 



dx / dw ^(^p^{x) — g4>{x)Y + t^^Qtix)^ 



(50) 



Once again, we assume that is very small compared to an infrared cut-off 
frequency fi. 

Define creation and annihilation operators by 

dx 



l-A 



{\!}^\4>{x) + m{x)) 



so that [d„,d]„] = Smm and then introduce the Bargmann-Fock variables a„(s) 
and a„(s). (For the mode n = 0, we arc proceeding as though the field had an 
infinitesimal IR regulating mass e << ^, so that 




but we will keep this replacement implicit for brevity.) Introduce as well the 
primed variables a'„{s), a'„{s) needed to express a mixed state of the field. 



* Attaching an independent bath to every point is an approximation, valid for 
sufficiently long wavelengths, as discussed in Section IV. Similar models, at zero 
temperature, have been discussed in References [12] and [13]. 



In this section we once again use a path 
to derive a propagator, similar to that in (14) 
scalar field: 

p[ai,a'J;t] = {a{,,| mWl) 

dald£,k d 



n 



27ri 



k,l= — oo 

X J[o:l,ai\ 

where \a'^) is shorthand for the infinite tens 
likewise for the bra. 

The propagator J may be computed frc 
to (17): 

J[ai,a';,ai,a'J;t] = Jv4,V, 

where the field variables 4>{x, s) and Il{x, s) a 
variables according to 

dx 



J -A 



a^(s) = 



(I 



dx 



-A2-\/|m|7r 

The primed version of (53) defines 4>' and 11' 
in (52) are as follows: 

dx -im^ fifn-n 



a„(o)= r 

J -A 



OCm{t) = 



dx 



dx 



aL{0)= [ - 

J-A2^\m\TT 

aL{t)= f- 



A 



.2^1^ A 
The boundary term B in (52) is most bi 

oo 

^=T- ^ (Qi,a„(t)-ha„(0)a: 
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where the boundary values of the a's may be found from the boundary values of 
the (f>'s and II's using (53). As in Section II, the action term in the path integral 
is of the usual pq — H form (up to a boundary term), minus a like term in primed 
variables: 

A=^ [ ds I da; (n</> - 011 - n'</>' + 0'n' - + n'^ - 9^^^ + 9„</>'^) 



-n+ 



n+n_ - 8^+84-^ , 

(56) 

whore wc introduce the linear combinations 0± = <f> ± <f) and Il± = 11 ± 11 . 

Because the heat baths coupled to the field at each point are independent, 
the influence phase V is simply an integral of influence phases similar to (21) 
(with I{lu) set equal to one and some simpliflcation via integration by parts, as 
in Reference [11]): 



ds ds' / dx(f>-{x,s)4>-{x,s') 
Jo J-L 

ds / dx^+{x,s)<p-{x,s) 

J-L 

2 

^ j da;0+(a;,O)</._(a;,O) . 



awwcoth— e ^ 



Trg_ 

4 



(57) 



The Euler-Lagrangc equations for this extrei 



2 
2 



4>-ix,s) 


= n_(.T,,s-) 


4>+{x.s) 


= n^(x.s) 


U-{x,s) 









„2 /■* 



ds 4>-{x, s') 



J 



Differentiating (59) and (60), and employin 
obtain the second order equations 



da;2' 
d^ 



ng 



<^+-^<P+ + ^-e{L-\x\)<p+ 



e{L-\x\)^- = 



2 



2 rt 
I C 




B. Solving the equatl 



Since we will once again be interested in the case where is small, we will 
neglect the last term in (57), on the grounds that it will produce only negligible 

perturbations on the solutions to the field equations. The black body is here 
assumed to have a uniform initial temperature and the limits on the integral 
over uj are meant to be taken to infinity only after all other integrations have been 
performed. 

The exponent of the integrand in (52) is thus quadratic, and the propagator 

J is proportional to the value of the integrand when this exponent is extremized, 
subject to the boundary conditions. Once again, A + V vanishes on shell, and the 
propagator is given by 



J[ak,ai ,ai„,a,l;t] = —e 



(58) 



We will find solutions to (63) and (64) tl 
and space. We begin with (63), and consider 



-{x,s) 



where A; is a complex frequency to be determ: 
satisfy 

2 ^ 4 

To be a saddlepoint of the path integi 
must be everywhere. We must therefor 
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approach x = ±L from either side. It is easy to see from the discontinuity in the 
wavenumber implied by (66) that there will be reflection at the boundary, with 

2 

reflection co-cfHcient of order This implies that there is really no such thing 
as a perfect absorber, since below some threshold frequency there will always be 
signiflcant reflection. Nevertheless, if we restrict our attention to field modes with 
frequencies above an infrared cut-off Q,, such that << Q,, then we can neglect 
reflection, and consider our black body to absorb perfectly at all frequencies of 
interest. Satisfying the periodic boundary conditions at a; = ±A forces k to 
assume a discrete set of values. Imposing these constraints on L±, M±, and R±, 
we find the most convenient parametrization of the general solution, above the 
IR cut-off, to be given approximately by 

<P-{x,s)= J2 e''^"(B„/„(cc)e-*l'^l^+S„/„(-cc)e*l'^l=)e^' + 0(|J), (67) 



where 



fn{x) = f{xsgn(n)) 
fix) 



( g-(r-HT)ig-7^ 



g(r-H7)ig-72: 



- A < a; < -i 

- L < X < L 
L < X < A 



(68) 



The attenuation co-efficients 7 = 2^ ^ and F = ^^-j^ are defined so that / is 
indeed periodic. 

The general solution to the homogeneous part of (64) is of the same form 
as 4>-{x,—s). The driven solutions, using the source 4>-{x,s') as given by (67), 
may be found by the usual Fourier technique. Once again discarding from this 
particular solution terms which are actually free solutions, we obtain the general 
approximate solution 

^+(a;,s)~ e*Tr='(s„/„(_a;)e-*lTf I" +B„/„(a;)e^l'rN)e-T^ 



- e''^"coth^^(A„/„(a;)e-'l'^l'' -A„/„(-x)e'lTl'')g7'' , 

\^\>n 

(69) 



where the "~" means that we have neglecte( 
all modes with frequency less than f2. 

Using (69) in (59) and (60), and droppi 



ill_(a:,s) 



E I 



\^\>n 

I |e^T-coth^(^„/„(x)e 



l¥l>n 



With regard to the neglect of modes below 

modes still have solutions of the form (65), 
non-negligible reflection. When we return to 
(53), therefore, their contribution to the Foi 
be negligible, for \ '^\ > Q,. Hence the iufra 
from the absorbed modes, and as long as we i 
the cut-off, we can use (69) and (70) without 
of the infrared contributions. 

This last argument hinges on the fact t 

with \ri\ « \k\ so that the growth or deca; 
over very many wavelengths, is primarily com 
close to k. This statement may be made moi 



d„ 



i(n—m)- 



. (n-m)7i 

' X — 



2iA 

where we take the upper (lower) of the ± sig 
quasi-diagonal — predominantly concentrate 



the diagonal. If we let 



while keeping F 



delta. In fact, the case where 75- e'^^ is not si 
the large factor e^^ will cancel in our flnal rei 
off-diagonal elements of dmn still suppressec 
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reflection terms ignored in (67) remain negligible, despite the fact that some will 
have amplitudes of order e^^ . 

The quasi-diagonaJity of dmn will be of crucial importance in the remainder 

of this paper. We will use it to construct an approximation whereby modes of 
all frequencies above the IR cut-off are slowly modulated, so that they effectively 
vanish in large regions of space, but retain well-defined wavelengths in the regions 
where they do not. It is in this manner that wc will identify thermal radiation 
of all frequencies propagating outwards from our model black body. The "no 

2 

reflection" condition that f^- << 1 thus turns out to have additional simplify- 
ing consequences far beyond making reflection negligible. As in Section II, it is 
weak coupling and slowly varying spectral density that lead to canonical thermal 
behaviour; the non-reflectivity of a black body is simply a sign that the weak 
coupling limit applies. 

We can now re-combine (f>± and II± into the Bargmann-Fock a variables, hav- 
ing effectively used the other variables to de-couple the Bargmann-Fock equations 
of motion. Combining (69) and (70) with the definition (53), assuming | ^| > fl, 
and then using |njiim„ ~ \m\dmn, we deduce the following approximate solutions 
for the Bargmann-Fock variables: 

f3\n\n 



amis) ■ 



n= — oo 

CO 

E 



-dmn(cOth ■ 



2A 



l)e'"An+dl„e-'"B, 



-d„„(coth + l)e^M„ -f- dl„e-^"B, 



rcoth -M)e^M„ d„ 

2A 



d1^„ (coth 



P\n\n 
2A 



l)e^M„ + dr. 



'Br 



(72) 

Once again, the dotted equality sign indicates that these equations are valid up 

2 

to corrections of order ^ . 



C. Boundary conditions 



We must now constrain An,A„,Bn, S„ to meet the boundary conditions of 



(54). We begin with A„ and B„. The const 
decoupled by defining 



An = 



coth 



1 y 

2A ^ ^ fe=-oc 

OO 



Using the quasi-diagonality of dmn, we can s 
mately by 

c.(i+c,f!))./;j 

since the matrix product dm.nd~l will, fo 
butions where m, n, and k all have the sam 
fn{x) — » fm{x) or fk{x) with negligible error 

The inverse of d^n is thus nearly diag( 

d*^^e~''"^ (since if /3 is large enough for Vf} 
and the whole term is negligible). This allows 

k.n 

a'£ ~ \/ \m\ir'^^^{dm.ne~^ "^^ 'd'^ e'*'* — d\ 



We can expand (75) into 



fc,n 



t J -A 



da^^i(m-k)^^i\k\ 
2A 
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where Gmix, t) = e'^ fm.{x)fm.{—x + t sgn{m) ) (continuing / periodically beyond 
[—A, A]). These last two lines are obtained via another quasi-diagonal approxi- 
mation: because the dominant contributions to (75) come from the terms where 
m, n, and k are all within FA of each other, and therefore (for | ^ | > Q) all have 
the same signs, we can switch /„ — > and \n\ — » nsgn(m) or nsgn(fc), and 

2 

incur errors of order ^j- at worst. This type of approximation will be used re- 
peatedly hereafter, drastically simplifying several matrix inversions. In effect, the 
only reason not to treat dmn. as a delta function is that, for large enough t, e'^* 
can vary significantly over its otherwise negligible width. For sufficiently small 

2 

jY, this quasi-diagonal approximation is accurate and natural. At the quantum 
level, the quasi-diagonal approximation is the ideal black body approximation. 

Applying this approximation again allows us to invert (76), obtaining 
_ a/Mtt J -a 



E 



(77) 



K 2A 



G,n(x,t) 

Note that the function Gm{x,t) = G(x sgn(m) ,t), where for t > 2L 



G{x,t) = f{x)f{-x + t)e 



For t < 2L, we have instead 



■yt 



G{x,t) = fix)f{-x + t)e 



It 



Ay* 



x G [-A, -L] 

X G [-L,L] 

X G [L,t- L] 

X e [t- L,t + L] 

a: G [t -I- L, A] 

X G [-A, -L] 
X G [-L,t- L] 
X £[t- L,L] 
xe[L,t + L] 
xe[t + L,A] 



(78) 



(79) 



As long as L and t are both sufficiently large, the diflferences between the two 
cases (78) and (79) turn out to be insignificant. 



Applying another quasi-diagonal approx 
(77), we obtain the boundary values 



, A — \ " I dx ^ 



\k\n y-A2A 



In a similar manner, we also obtain 

, A — Ofc I dx'^ 



|A;|7r i-A2A 



at, I ax 



V^A^i-A' 



2A 



Using (80) and (81) in (55) and (58), we 
J to be given by 

i -li -f If , 



^ exp^(af aLT<„ + a'ia'lTi 

l,m 
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where 



TL 



^ 2A 



2A 



G,n{x,t) 



L 

i: 



(83) 



2A 
2A 



These matrices are obtained using one new kind of quasi-diagonal approximation, 
which may be illustrated by writing 



Tin 



G^{x,t) - ^^^(^ t) 



A 2A 



,|m|„ 



(84) 



.1^ 



Gm{x,t) 



2A 



,\1V 



t) 



by substituting x ^ x — (x — t sgn{m) ) and observing that Gm (—x + t sgn{m) ) = 
Gm{x,t). (The last step in (84) invokes the kind of quasi-diagonal approximation 
already introduced in (74).) 

Equations (82) and (83) provide us with the propagator for the density 
matrix of the modes above Cl of the masslcss scalar field, in the presence of a 
black body at inverse temperature /3, initially uncorrelated with the field. The 
integrals in (83) do not seem to have closed form evaluations, but since T, T' , A, 
and E are all matrices formed by adding and multiplying matrices which axe quasi- 
diagonal, they must themselves also be quasi-diagonal. (The inverse of a quasi- 
diagonal matrix is quasi-diagonal; one way of showing that this statement applies 
in our case is to consider the massless field as the continuum limit of a lattice, 
in which case all the matrices become finite dimensional, and the proposition is 



then elementary. One might be concerned th 
diagonal matrices could cancel when they ai 
of the matrices in (83) can be computed, a 
quasi-diagonality of the matrices in (83) wi 
body propagator J, by calculating some sigr 
9'. 



C. Physical interpretation 



For convenience in explaining the signi 

to some important spatial regions; these regi 
overlapping intervals [—{L + t),L] and [—L 
left and right shadows of the black body. W 
be the smallest length for which e 2 » is n 
[-(L + t),-(L + t-\)] and [{L + t - A), (L 
while [-{L + t- X),-L] and [L,{L + t-X)] cc 
[—L, — {L — A)] and [{L — X),L] are referred tc 
(L — A), (L — A)] forms the core. As long 
the differences between the two forms of G{x 
no effect on our identification of these region 
insignificant in (83). 

With these relevant regions in mind, a 
gator is immediately suggested by the form 
observe that Gk{x,t) — vanishes for 

body; it has support within the left shado' 
shadow for A; > 0. therefore seems to 

firom the black body, and Amk seems relatec 

into it. As well, (1 — e '^~^)[Gk{x,t) — ^^j-jj 



umbral and black body zones, but equals unil 
are therefore suggestive of ordinary, free field 
disconnected from the black body. 

These impressions can be confirmed bi 
involving J. Firstly, we can consider the 
evolves under J, and show that all iV-point 
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are the same as those for uniform, unidirectional thermal radiation, as long as all 
the N points are within the umbra of the same side of the black body's shadow. 
Secondly, wo can show that the uniform thermal state at the temperature of the 
black body evolves under J' into itself. Finally, we can set the temperature of 
the black body to zero, and obtain the "quantum wall effect" whereby quanta are 
absorbed by the black body, without reflection. These three special cases suffice 
to illustrate that (82) is indeed the correct quantum mechanical description of 
the effect of a black body on a massless field. 



The final state evolving from initial vacuum is found simply by setting 



=0 in (82): 



(85) 



We derive this by using quasi-diagonality to 



The expectation value (87) can be cor 

states of uniform, isotropic thermal radiati( 
radiation travelling in one direction only. Thi 
to that of (88), which are 



= E 



- 1 



The correlation functions for products of operators 4>{xi)...4>{xN)Il{yi)--.tl{yM) 
may be formed from the expectation values 



HE' 

n—l \ m 



2A 



HE 




(6) 



(86) 



/ n— 1 \ m 

where Umn,Vmn are the factors (■'^x^)*' applicable to either or H. This ex- 
pectation value may be computed by twice completing "squares" ^ in a Gaussian 
integral: 



n 




I\ IVl 



K„=K„=0 



where 



E- 



VkmUkn 



(87) 



(88) 



l-e" 



t Actually, bilinears in Bargmann-Fock variables are not squares, since a and a 
are not complex conjugates, but the procedure is exactly analogous to the usual 
method of manipulating quadratic exponents. 



-(j/m-^n) 



'-{ym- 



fc=o e- A - 1 

By comparing (88) with (90) and referring t( 
far as measurements made outside the shade 
the state represented by po [ai , Um ; t] is indis' 
measurements made in the left (right) umb 
state where only left-moving (right-moving) ] 
for measurements inside the core of the body, 
thermal radiation. (The asymmetry in (87) 
spurious, since Gi {x, t) differs appreciably frc 
the thermal kernel defined by the sum over 
j^hy)- The asymmetry therefore contribute: 
disregarded to be consistent in our ideal blac 

For our second illustration, we compute 
the thermal state at the black body temper; 
operator 

1 V~ 

p,3Km,^m;0] = ^exp2_ 

m 

and using (51), integrating by completing th 
pp[ai,am;t] = ^ exp^Q;i[T • (e'^~ 
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We can approximate 



,|m|„ 



i(i-m)^ 



G,n{x,t) 



Gm{x,t) 



G™(cc,t)(e^ 



1) 



(93) 



which leads to 

^ ^i(m-i) S! 



/.A 

;n(m) ^ / 



2A 



-(1-e- 



Gi(x,t)-e- 



,|m|„ 



(94) 



/da 
.a2A 



■(1-e- 



,|m|„ 



) 



We make this last step by again changing variables x — + — (x— t sgu(m) ) (wrapping 
around periodically in the domain [—A, A]). It immediately follows that 



which establishes the fact that 



(95) 



as far as modes above the IR cut-off Q. are concerned. The canonical ensemble at 
the temperature of the ideal black body is indeed the equilibrium state. 

Finally, we can consider the evolution of a one-particle initial pure state, 
in the non-emitting case where /3 — » oo, so that 0. Changing to the N- 

particle basis by extracting the co-efficients of a^a^ and cJ^ol^ in »7|/3_,oo, we 
find that the final state evolving from the initial state 

E, i , /* I ^TT \ / mTT I 



where j^) denotes the one particle state with momentum has the density 

operator 



Pi 



=m(oi+(i-F)5:^/v.i-i^)(^i 



(96) 



where |0) denotes the vacuum and 



In the limit /? - 



oo, we have 

Tim „ 



This ensures that an initial state localized wit 

of the black body evolves into a state in whic 
above the IR cut-off are negligibly excited. A 
black body evolves into one in which no mod 
excited. This then implies that the probabi 
black body is negligible. Note that, when (3 - 
This shows that the total probability is co: 
particle being present decreases from one. 

Our model and the quasi-diagonal appi 
it therefore exhibit absorption, thermodynar 
diation. We have presented, in an idealize 
a description from quantum mechanical firs 
macroscopic black body with a massless qu 
closely resemble that of light. 



IV. DISCUSf 



A. Summa 



The use of influence functionals to desi 
unobserved environment has not received its 
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theory. Among second order effects, there can be some that are enhanced by large 
parameters that do not appear at lower orders. These are very important, since 
they arc apparently responsible for the significant phenomena of quantum mea- 
surement and thermal dissipation. As discussed by Feynman and Vernon in their 
classic paper introducing influence functionals[9], independent harmonic oscilla- 
tors can be used to describe any environment, as long as we are only interested 
in effects up to second order in the coupling between environment and observed 
system. To this order, a general quantum mechanical environment having allowed 
transitions with energy differences ujk is equivalent to a set of harmonic oscillators 
with those natural frequencies, with effective coupling strengths proportional to 
the matrix elements of the interaction Hamiltonian^ [9] . 

This result is derived in the Appendix. To first order in the coupling, a 
generic quantum environment may be treated as a classical source; the first 
improvement beyond this treatment, valid to second order in the coupling, is 

to model a generic environment as a collection of independent oscillators. Our 
treatment of matter as a collection of baths of harmonic oscillators is therefore 
not merely a gross idealization, but is actually accurate up to second order in 
perturbation theory. 

Prom another viewpoint, though, our model of a free quantum field interact- 
ing with first-quantized matter can actually be considered as a way of approach- 
ing some non-perturbative physics. The atoms and molecules represented by heat 
baths are in reality bound states involving the electromagnetic field. Interacting 
fields can effectively have many more degrees of freedom than free fields, because 
their Hilbert spaces include arbitrary numbers of bound states. These bound 
states couple to the unbound modes, and can have significant effects if they are 
present in macroscopic quantities — as in a black body of macroscopic size. 

That these effects can be significant despite weak coupling can be starkly 

demonstrated by calculating the energy radiated by a black body. From our one 
dimensional model in Section III. we find (in the limit where the infrared cutoff 



One might be concerned that, while the dimensionful set the absorption scale 
A in our one- dimensional calculation, in actual spacetime the electromagnetic fine 
structure constant a^m is dimensionless. In fact, though, the effective couplings 
which appear when realistic matter is approximated by harmonic oscillators must 
be proportional to dipole moments of molecules, etc., and so will still have the 
correct dimensionality. 



is removed) that the outgoing radiation carri 

6/32 &n 

making % and the Boltzmann constant ks & 
ature. This is the one-dimensional, spin ze: 
law. 

The Stefan- Boltzmann constant a is indepenc 

yet it governs electromagnetic radiation gene 
of Section III we can now understand this 
radiation occurs in situations with microscop: 
one must perturb separately in ^5- and e~' 
quantities are small is qualitatively different f 
to zero. The Stefan-Boltzmann constant mi 
"invisible" charge-dependent factor of a for 
vanish if the electromagnetic coupling were 
negligibly different from unity for macroscop 

The condition that the black body has n 
spatial growth and decay of the field solutio 
their oscillatory frequency. This in turn impl 
approximation, whereby modes of a given free 
of the same frequency, but with spatial modr 
large enough distances. An approximation of 
whenever we speak of light of definite frequei 
finite regions of space. The prevalence in na 
that justify such concepts is ultimately due to 
The vanishing refiectivity of a black body 
indicator of the underlying physics of weak c 

Compared to the physics of weak coup] 
der effects by macroscopically large paramet 

The effective coupling may be even weaker, i: 
in deriving the independent oscillator heat 
theory. The density of degrees of freedom ir 
does not seem unreasonable to conclude tha 
way, as medium-dependent as classical optics 
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familiar physical phenomenon. Yet it also permits a powerful type of approxi- 
mation, elucidating "qualitative" properties which, like black body radiation and 
quantum measurement, have traditionally been treated by axiom rather than by 
perturbation. The opacity of a black body can be considered as an indicator that 
a "large parameter" approximation is appropriate. 

With two very different, very drastic approximations applying at once, it is 
little wonder that a non-rcflcctivc, opaque object is such a clean theoretical sub- 
ject, and why its behaviour is so universal (depending only on temperature). The 
highly non-trivial implications of total absorptivity, first deduced from thermo- 
dynamic axioms, have now been traced to their microscopic origins in quantum 
theory. The results of our seemingly very idealized model, with uncountably 
many independent harmonic oscillators coupled to the field at every point, are 
thus quite universal. Is this model then more realistic than it at first seems? 

B. Justifying the model 

Apart from the reduction to one dimension, this model is in fact quite real- 
istic. In a disordered collection of many molecules there are very many allowed 
transition energies, and the continuous spectrum /(w) becomes a reasonable ap- 
proximation to a discrete spectrum of equivalent oscillators, that will not break 
down until one considers a time scale on the order of the inverse of the spectral 
spacing. Since a continuous spectrum of harmonic oscillators has infinite heat ca- 
pacity, one would expect this time scale to be associated with heating and cooling 
of the black body. These important phenomena are not described by our model, 
because of our assumption of a continuous spectrum. 

In Section III, however, we assumed that there is an independent heat bath 
interacting with the field at every point in the black body, so that our model is 
continuously dense in space, as well as in frequency. This is a valid approximation 
for a discrete model, as long as we consider wavelengths of light long compared 
with the lattice spacing. We axe therefore certainly only modeling the effect 
of matter on light with wavelengths above the X-ray range. In fact, we are 
assuming that each heat bath represents a grain or clump of molecules or atoms 
that effectively couples to the field at a single point, and so must be of a size less 
than the wavelengths of interest. For wavelengths as small as the near ultraviolet 
range, this permits us to have on the order of a thousand atoms per heat bath. 



This number of degrees of freedom must sti 
bath cooling time scale much longer than ou: 
a conflict between our assumptions of contin 

We can estimate the cooling time scale fc 
its classical thermal energy ^NksT by the ra 
sical black body of radius two to three nano 
law (100), we find the cooling time to be 

ISOOfeBT ^ . 

Tr- ^ ~ 10 

We can bound phenomenologically by assu 
infrared, at around 10^^ Hz, and by knowing t 
a thickness of around 10"'' metres, giving a ra 
seconds. For temperatures ranging into the 
have Tc >> g~^- The continuous spectral der 
not break down for our local heat baths until 1 
are well established. At this point issues aris 
the medium, that are beyond the scope of 
relate to a further limitation on our model: 
independent of all the others, we are explici 
effects involving longer range correlations in ; 
weak dielectrics, on the other hand, should 
model. 

The implicit infrared and ultraviolet cu 
in confiict with the ohmic spectral density 
stancy of the spectral density is not essential 
slowly varying leads to similar results 

constant r(tj). As long as e"^^^'"^^ is negligil 
absorptive. 

As far as infrared to near ultraviolet ligh 
presented in this paper is actually a fairly a' 
composed of closely packed grains, whose si; 
order of a few nanometers, each of which cons 
fine-grained deposits of carbon atoms occur 
used for making actual black bodies black 
makes a rough surface — a feature we have 
turns out to be about soot. 
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C. Issues raised by the model 

The problem we have studied presents a new face of the basic quantum 
dilemma regarding localization in position and momentum. One may often think 
of a black body as thermally exciting outgoing field modes, and absorbing incident 
modes, but this is misleading. The positional and directional localization involved 
in identifying outgoing and incoming radiation is more subtle than the one-to-one 
mapping by which we associate Fourier modes with particle momenta. 

This subtlety can be illustrated by sketching the extension of our work to a 
black body which is a solid sphere in three dimensions. Choosing spherical polar 
co-ordinates, the radial eigenfunctions of the free field Hamiltonian are spherical 
Bessel and Neumann functions. While these can be combined into spherical Han- 
kel functions, describing incoming and outgoing waves, the Neumann functions 
diverge at the origin, and the only physical modes are pure Bessels, represent- 
ing standing waves. There are not enough field degrees of freedom, then, to be 
decomposed into orthogonal outgoing and incoming modes. 

The black body environment surrounding the origin changes the path in- 
tegral saddlepoint, and modifies the Euler-Lagrange equations. Taking the zero 
angular momentum mode as an example, we find /o(f) = ^-—^ — - — " ' — - ap- 
pearing in solutions, instead of 2HL*i, Radial functions that arc predominantly 
incoming or outgoing (outside the black body) will therefore appear in the three- 
dimensional analogue of (83). But the number of orthogonal field modes is not 
suddenly doubled, and each of the incoming and outgoing sets of eigenfunctions 
could in fact be expanded in terms of the other. There is only a half-line k £ [0, oo] 
worth of radial degrees of freedom to be described. 

This point is one more illustration among many that the wave-particle du- 
ality of quantum fields is not completely explained by referring to the expansion 
of the field operator in plane waves. While it is convenient for scattering pro- 
cesses in negligible gravity, this notion of particles fails to be physical in curved 
spaces or accelerating frames; it also requires elaboration (via our quasi-diagonal 
approximation) to describe emission and absorption from macroscopic bodies at 
rest in fiat space. A really satisfactory formulation of particle states would be of 
great value, and it has yet to be found. 

A second issue is raised by our zero-temperature, purely absorptive limit. 
The "quantum wall effect" blocks field excitations from passing through the black 



body, but it does not affect the ground state 
boundary condition one typically assigns to a 
there is no Casimir effect, at zeroth order ir 
of non-conducting dielectric[13]. The case 
this paper, should also admit a similar kind 
point heat baths replaced by a gas of free 
would be interesting to study such a model 
effect without simply imposing boundary con 
important phenomenon of thermal radiation 

Finally, we could relax the assumption ( 

implied infinite heat capacity and discarded i 
cooled to low temperatures, we would expect 1 
from the field would be re-emitted, as it appr( 
problem would also provide an interesting coi 
of a black hole through the Hawking effect, 
to undergo radiative heating rather than co 
trapped information is common to both blac 
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APPENDIX: INFLUENCE FUNCTIONAL FOR 
A WEAKLY COUPLED ENVIRONMENT 



This path integral may be computed e 
HE{p,q) = +a;2g') and A{p,q) = q. W( 



Consider an environment represented by the single degree of freedom g, with 
conjugate momentum p. Let this environment be coupled to an observed system, 
such that the total Hamiltonian operator is 

H = Hs + HE + eA{p,q)Q. (102) 

The coupling constant e is assumed to be small and dimensionless. Hs and Q 
are operators on the observed system, while He and A{p, q) are operators in the 
unobserved sector. 

Let 17^) denote cigenstatcs of the environmental Hamiltonian, such that 
He\K) = Ek\K), assuming no special form for He{p, q). We can then in principle 
define the q and p representation wave functions tpKiq) = {q\K) and ^k(p) = 
{p\K}. Let the mixed initial state of the environment at time s = be described 
by the density operator R =:R{p, q):, where the colonss imply normal ordering by 
placing all p's to the left of all q's, so that R{p, q) is the Wigner function {p\R\q). 
R has the matrix elements 

{K\R\L) = Rkl 

in the basis of energy eigenstates. The Wigner function may be expressed in 
terms of these matrix elements using wave functions: 

Rip, 9) = *K(p)J^x4o)V'a?) ■ (103) 

K,L 

For the environment with this Hamiltonian and initial state, the influence 
phase V[Q, Q'\ is given by the path integral 

giV[Q,Q'] ^ j dpfdq'fdq,dp[R{q.,p[)e'^''f'''f 
•'n;0 Jp'fi 

(104) 

Here A' and H'^ are short for A{p' ,q') and HB{p',q'). 



VIQ, Q'] = ^J^ dsj^ ds' (Q - Q')s 

X ((<? - <?')»' coth cosuj{i 

where s and s' are time variables, (Q — Q'), 
time. 

The influence phases for any number of ' 
environmental degrees of freedom simply ad 
systems are in a thermal state, each will be s 
related with the other systems. Consider, the 
a number N of independent harmonic oscilla 
ral frequencies u)k and coupling strengths ek- 
thermal state and uncorrelated with the obsei 
functional of the form 

/•oo rt PS 

y[Q,Q'] = j2 / d^Iis^) ds ds'{Q- 
" Jo Jo Jo 

X ^coth ^(Q — <5')f.' cos w(s 

where 

JV 

fc=l 

As we have seen in Section H, above, 
influence phase of the same form as (106). 

Returning to the environment with onlj 
consider the case of general i?E(p, g) and A 
small enough, equation (104) may be approxi 
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in e, we have 



■'^'l = j dpfdq'fdq,dp[R{q,,p'^e-i'''''f 



X / VpVq I VpVq e'^ 
p';0 



JJds {(p<i-p'q)-(HE-H'^)) 



X ^1 - j ds{AQ- A'Q') + 0{t'] 
= 1-| / ds^RLMAMLe-''''^^%Q-Q'), + 0{^) 

•'0 T. M 



(108) 



: exp 



(-'-lj\{A{s)){Q-Q%)+0{^) 



In (108) we use the following notation: 

Ar^M = {L\A\M) 
El — Em 



n 



{A{s)) = RlmAmlb 



When we recall that the influence phase appears in any path integral over 
observed degrees of freedom as an additional term in the exponent, 

g^(s[Q]-S[Q']) _^ gt(siQ]-S[Q'] + V[Q,Q']) ^ 

we can see that the first order term in (108) clearly constitutes a perturbation 
to the observed sector action S[Q], with the expectation value {A{s)} acting as a 
classical external source. 

Since the first order term leads to effects that are already well studied in the 
regime of closed system quantum mechanics, we generally wish to concentrate 
on the specifically open quantum behaviour produced by the higher order terms. 
If we restrict our attention to environments that are approximately in equilib- 
rium, so that Rmn = ■^e~^^'^ 5mn, and assume that the expectation value of A 



X / VpVq / Vp'Vq'e' 
9«;u ''p'i'O 



vanishes^, then the first non-trivial term in (1 
giv[Q,Q'l ^ j dpsdqjdqidPiR{q,,Pi)e-^''i' 
j-Pf 

[ ds {AQ 
\ [Jo 

•'° L,M 

n ^ Jo Jo r,^,^ 

X ((Q-g') 



' coth 



-iiQ + Q')s'S 

This second order term cannot be repr 
action; it contains qualitatively now behavioi; 
system. Using the smallness of , we can re- 

VIQ,Q'] = ^ J dwG{u,;l3)J^ dsj^ ds' [q{s 

X ^coth ^[g(s') — q'(s')] cos w(s — 

where the effective spectral density is 

G{uj;l3) = ^ <5('^ - ujlm)Almj 

L>M 

The restriction of the environment to 
ally irrelevant to the above derivation, and n 

t By considering a slightly generalized therme 
order external sources together with the se 
paper. 
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convenience. Comparing (110) and (106), therefore, we see that any weakly cou- 
pled environment that is initially uncorrelated with the system and in a thermal 
state is equivalent, at a fixed temperature, to a bath composed of independent 
harmonic oscillators. It is worth emphasizing that the effective spectral den- 
sity G(a»; /3) can be temperature dependent, as this fact seems to have escaped 
comment heretofore. 

In fact, (110) can also apply in a situation where the weak coupling ap- 
proximation is not £is straightforward, namely the case of a large number A^" of 
independent {i.e., uncoupled and initially uncorrelated) environmental degrees 
of freedom, each of which is weakly coupled to the observed system, but not so 
weakly that the total effect of the environment has to be small. In this case the 
Hamiltonian may be written as 



H = Hs + 



N 

E 

fe=i 



{Hk + ekAk{pk,qk)Q) , 



(112) 



where all of the are much less than unity. 



Assuming again the initial thermal state, with vanishing expectation values 
{Ak{s)), we find 



JVIQ,Q 



'1 = J dpfdq'fdqMe-i'''f'fR{qi,p'0 
I VpDq I VpVq 



JV 

X I I e 

fc=i 



ds (AQ - KQ') 



(113) 



ds I ds' I dujGk{oj;P) 
Jo Jo 



k=l 

X ((Q - Q')s' coth cosa;(s - s') 

- i{Q + Q')s' sma;(s - s')) + 0(£?)) 



Here Gk{ijj;f3) is the effective spectral density for the fcth environmental degree 
of freedom, defined similarly to G(w;/3) in (111). 



In the case of large N, we can now now 
1 /N. The leading order term may be found I 
taking the limit iV — » oo. In this limit, the p 



lim r\i'-^4h 



and so we see that the influence phase for tl 
of weakly coupled degrees of freedom is givei 

V[Q,Q'] = {l + 0{ek) + Oij^)) X ^ J 
X ((Q - Q')=' coth ^ 

- KQ + Q')s 

where the effective spectral density is 

JV 

Equation (114) is, to leading order in tk an< 
(110) and (106). 

For large enough N , quantities such as 
derivations of wave function collapse [6] — m 
the Efe are all small. It is important to no1 
coupling of the observed system to each em 
be weak, there may be so many such dcgr 
the system can be large. Equation (114) s 
generic environment by a bath of independt 
in this case. We have therefore confirmed tl 
oscillators provides an adequate model for a 
large environment, capable of inducing drasti 
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